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Abstract 

We examine the Cerenkov-like emission of e^e~ from muon super- 
luminal muon neutrinos assuming a quadratic energy dependence of 
the neutrino velocity arising from Lorentz violating interactions. We 
find that with the OPERA result for the neutrino-photon velocity dif- 
ference , the decay length for the process i/^ v^e^e~ is 17,039 
km which is much larger than the OPERA neutrinos path length 
of 730 km. We also calculate the pion rate for super-luminal out- 
going neutrinos, and we find that the deviation of the pion decay 
length from the standard Lorentz conserving case at the OPERA neu- 
trino energy is 2%. We conclude that if the muon-neutrino velocity 
has a quadratic energy dependence, then OPERA result is consistent 
with non-observation of forbidden neutrino decays and large devia- 
tions from the standard pion decay lifetime. 



1 Introduction 



The OPERA experiment has recently claimed to have observed neutrinos 
traveling faster than light [1]. A possible explanation for such superluminal 
neutrinos comes from Lorentz violating interactions. An important phe- 
nomenological constraint on super-luminal neutrinos is from the observation 
of Cohen and Glashow (CG) [2] that super-luminal muon neutrinos will lose 
energy via a Cerenkov like emission of an e^e~ pair. 

The OPERA result find that muon neutrinos of average energy 17.5 GeV 
traverse a distance of 730 km from CERN to Gran Sasso with a velocity 
which exceeds the photon velocity, 

5{E = 17.5GeV) = ^'"'^ ~ = (2.48 ± 0.28{stat) ± 0.30{syst)) x 10"^ (1) 

This is comparable with an earlier measurement of muon neutrino velocity 
by MINOS [3] who found that muon neutrinos of average energy 3 GeV 
traversing a distance 730 km exceed c by an amount, 6{E = 3GeV) = (5.1 ± 
2.9) X 10^^. This is in contrast to the neutrino observations from supernova 
SN 1987a [4, 5, 6] where over a flight path of 51 kpc, the neutrinos with 
energy in the band (7.5 — 39) MeV all arrived within a time span of 12.4 sec 
and the optical signal arrived after 4 hours of the neutrino signal (consistent 
with prediction of supernova models) from which it is inferred that 6{E = 
15Me\^) < 10~^. This implies that the to be consistent with all observations, 
the neutrino velocity is energy dependent [7, 8, 9]. 

Horava-Lifshitz theories [10, 11, 12] provide a framework where theo- 
ries are made renormalizable by the introduction of Lorentz violating higher 
derivative terms in the Lagrangian. Models of Lorentz violation which can 
give an energy dependent neutrino velocity are discussed in [13, 14]. The 
Lagrangian for the Lorentz violating neutrinos is given by 

C = ^(^zl/)-m-^{u-Dr-'-^^{u. Dfiu ■ 7)) Pl^ (2) 

where m*^ is a fixed four-vector which represents a preferred frame, thereby 
explicitly breaking Lorentz invariance and the scale of Lorentz violation is 
determined by a large mass M and the dimensionless parameters ai and 02- 
Starting with the Lagrangian (2) the dispersion relation for neutrinos can be 
derived of the form, 

4 

E^=p^ + m' + 7^'p^ + '^ (3) 
where 77' = mai/M and 77 = 2a2- 
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The neutrino velocity from eq.(3) is given by 

dE 1]' Srip^ 

S = 1 ~ — H — 

dp 2 2 M2 



(4) 



At energies p ^ y/mM, we have rjp'^/M ^ rj' and the neutrino velocities 
increase quadratically with energy. 

It has been pointed out, the modified dispersion relations for neutrinos 
make the process u — )■ uff kinematically possible [2, 13, 15, 16] and the 
pion decay lifetime can get a sizable modification from the phase space of 
the outgoing super-luminal neutrino [17, 18]. 

We will assume that (a) energy and momenta are conserved in all Lorentz 
frames and (b) the energy-momentum relation for neutrinos is of the form 
(3) in the lab frame. For the electrons and other particles are assumed to be 
of the standard Lorentz invariant form Ef = pj + rn^- 

In this paper we compute the processes z/^ — )■ z/^e"'"e~ and vr — )■ /ii/^ 
assuming the muon-neutrino dispersion relation (3) in the two limits: 

1. E'^ = p^il + 7]') ( 5 is independent of energy and we ignore neutrino 
masses) and 

2. E'^ = p^ + T]p'^/M'^ ( 6 is quadratic in energy). 

Our results are as follows. Assuming E"^ = p'^{l + rj') we find that the 
lifetime of the z/^ — )■ t'^e"'"e~ is r = 882.9 km for E = 17.5GeV which means 
that more than half of the OPERA neutrinos should decay in a 730 km fiight 
length. We also find that the pion in fiight decay width decreases by 24%. 
These are large effects which are not observed and this rules out non-zero rj' 
as the source of the super-luminality of neutrinos observed at OPERA. 

On the other hand assuming E"^ = rip'^/M^ we find that the lifetime of 
the — )■ z/^e"'"e~ is r = 17038.6 km for E = VJ .bGeV which implies that the 
number of neutrinos is depleted by only 4.2% in the course of the CERN-Gran 
Sasso fiight. The decrease in the pion in fiight decay width is not observable. 



2 Neutrino energy loss by electron-positron 
pair emission 

2.1 Energy dependent neutrino velocity 

We assume the neutrino dispersion relation E'^ = p^ + rjp'^/M'^ in the lab 
frame and calculate the decay width of the process z/(p) — )■ v{p')e^{k)e^ {k'). 
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The amplitude squared for this process is given by 



|M|2 =32^1 [{p-k'){p ■ k) (1 -4sin2^vi/ + 8sin^^vK)] . (5) 
The decay rate of the neutrino is in general given by 



' (27r)32E^ (2vr)32E^ {2Txf2Ee 2E^ ^ ' ^ ' ^ ' 



Using 



d^k5{P)e{ko) 



we have 



1 f d^p' d^k' |M|2 , 

' 5{{p-p-k'f) (7) 



8(27r)5 7 El E'^ E, 

where we have performed the k integral and imposed the 5-function condition 
k = p — p' — k'. Without loss of generality we can choose 

p=(E,,0,0,|p|) 
p' = {El, Ip'I sin6',0, |p'| cos^^) 
k' = El{l, cos0sin6'i, sin sin 6^1, cos6'i). 

The argument of the S function in eq.(7) can be rewritten using the above 
definitions as follows 

ip-p- k'f = (^(IpI' - Ipf )(|P| - |P I) - IpIIp'I^') - DEI (8) 
where 

^ = ^(IpI' - Ipf ) + (IpI - W\)ol - \v\o^ + '2W\oe, cos<p (9) 

Here we have assumed that since we are dealing with high energy processes 
the angle of scattering is typically very small and of the order of rjp'^/M'^, 
thus dropping higher orders of ^, 6i and rj/M"^. From here on we shall use 
the notation p and p' to denote |p| and |p'|, the magnitudes of the initial 
and final state neutrinos respectively. Now we can rewrite the (5-function in 
eq.(7) as 

y (^^ " (m^^^' ~ ^"^^^ " ^'^ " ^^'^') 
Using eqs.(lO) in eq.(7) we get 

^ j p'dp' j de^ j E'jEi j del j d<j) 



5127r 
6 



{e:-{^,{p'-p"){p-p')-Pp'9^)d-')\^ (11 
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choice of momenta and remembering that k = p—p' — k', 

2 

W =8GI{E:e,p') {p - E'J9'9l + ^ { (P - K)0' + 2E'J9, cos + 9^ {p 



The |Mp, using our 
becomes 



/2 /3 /2 

p p'^ p'^ 

Now to fix the hmits of the 9"^ 



(1 — 4 sin 

^v^ + 8sin^^iy) (12) 

o iix me nnuLs oi me 9"^ and 9\ integrals we need to find their 
maximum values from the 5-function condition i.e 

DE'^ = ^^{p'-p"){p-p')-pp'9^ 

For the maximum value of 9 we set E'g = in the above equation 
have 

o2 V {P^ - P'^){P - P') 

'""^ M2 pp' ^ ^ 

And similarily setting p' = and the electron energy at its maximum i.e 
E'^ = p/2 in the 5-function condition we have, 



(13) 
so that 



r]p 

We make the following change of variables to pull out the factors oirj/M 
and p from the integrand : 



(15) 



M2 



(16) 



Using the above definitions in eq.(12) and substituting in eq.(ll) gives the 
rate of electron-positron pair emission as 



dx 



d9 



d9i 



(17) 



where 



(x- {1 + 9)x^ - + x^y 



9i{l-xf (^9i+x + 9ix 



f=-. Z ^= 

4 il + 9i{l-x) -9x -x^ + 2V99ixcos 

+2(1 + 9i)x^ + 2x^ + x*^ +9 (^9i{l + 9i) +x- 9ix - (1 - 9{) 
+2y^ (l - (1 - Oi - 09i)x - (1 - 9)x^ + (1 - 29i)x^^ cos 
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After numerically solving the above integral we get the following expression 
for the rate of electron-positron pair emission 



The decay width written in terms of 5 ~ (3/2)(?7p^/M^), is 



r = ^ 4 (1 - 4 sin' + 8 sin^ dw) (19) 



For OPERA neutrinos with energy E = p = 17.5GeV, the decay time is 



which means that the neutrino number reduces to the fraction N/Nq = 
exp(— 730/17038.6) = 0.958. A 4.2% reduction in the number of muon- 
neutrinos in the course of the CERN to Gran Sasso flight may be compatible 
with observations of the same neutrino beam by ICARUS [19]. 

At higher neutrino energies ~ SOOGeV and above, a neutrino produced 
in a collider will decay within the detector into hadrons and charged leptons, 



so it may be possible to test this dispersion relation at the LHC [20]. 

2.2 Energy independent neutrino velocity 

We calculate the rate for the process z/(p) u{j)')e^{k)e^{k') assuming the 
dispersion relation E"^ = + rj'p^ /M"^ which leads to a energy independent 
5 = 1]' /2 which is the same assumption as made by Cohen and Glashow [2]. 
And the expression for |Mp in eq.(12) now becomes 



r = - = 17038.6 km/c 



(20) 




(21) 



|M|2 =8GUE',E,p') [{p - E'MO' + 0^) + Iri {{p - pi 
+E'^eei cos 0}] (1 - 4 sin^ Ow + S sin^ Ow) 



K)ol 



(22) 



The 5 function in eq.(7) now becomes 



D'-^5 [e'^ - {r]\p - p'f - pp'e^)D 



) 



(23) 



where 



D' = 7]'{p- p) + pel - p{f + e\) + 2pieex cos 



(24) 
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Using the condition imposed by the 5-function we once again derive the hmits 
of the two angular integrals. Once more we put E'^ = and p' = in the 
5-function condition to obtain maximum values of 6 and 6i respectively. 

(25) 



PP' 

'ijmax 



- r] (26) 



And finally we change the variables of integration as before from p', 9 and di 
to X, 9 and 9i respectively with 

p' ^ xp, 9'^ ^ r]'9, 9l^r]'9i. (27) 

In this case the rate of electron and positron emission from a neutrino 
decay becomes 

r = ^ ^ V V (1 - 4 sin^ 9w + 8 sin^ 9w) ■ (28) 
Expressing the decay width in terms of 5 = rj' /2 we obtain 

r = ^ ^ Sy (1 - 4 sin^ 9w + 8 sin^ 9^) (29) 
For OPERA neutrinos with energy E = p = 17.5GeV, the decay time is 

r = Y = 882.9 km/c (30) 

which means that the neutrino number reduces to the fraction N/Nq = 
exp(— 730/882.9) = 0.437. A 56% reduction in the number of muon-neutrinos 
in the CERN to Gran Sasso flight can safely be ruled out [19]. This implies 
that the dispersion relation E"^ = p^(l + rj') to describe energy independent 
super-luminal neutrino velocities can be ruled out as pointed out in [2]. 



3 Pion decay lifetime 

3.1 Energy dependent neutrino velocity 

We calculate the pion decay width in the lab frame with a super-luminal 
neutrino in the final state. We assume the dispersion relation E'^ = {p^ + 
rjp'^/M^) in the lab frame. The amplitude squared for the process vr~(g) — )■ 
IJ'~{p)^n{k) is, 



m„ + 



7]k 



ml 



m. 



+ 2 



(31) 
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The decay width is then given by 



Ef, k 



+ 2 



S^(^q- p- k)S{E^ -E^- E,) [ml 



m, 



m. 



(32) 



Performing the d^p integral to remove the 3-momentum 5-function and 
writing = — k\^ + rnl the decay rate then becomes 



SvrK 



kdkd cos 9 
\J\q — + ^ 



-.5{E^ + \^\q-k 



ml - E^) 



2 , Vk^ 



m 



m. 



(33) 



Writing \q — k\^ = k"^ + — 2kqcos6, 9 being the angle between k and 
q, and E^ = k + rjk^ / {2M'^) we see from the argument of the 5-function in 
eq.(33) 

rjk^ ^ rjk^^ 



cos 9 



ml-ml + 2E^k + - 



{2kqr' 



(34) 



while the derivative of the argument of (5-function with respect to cos 9 yields 



d 



iE, + J\q-k\^ + ml-E, 



dcos 9 

Substituting this in eq.(33) we get 
, _ GlP^ml f dk 



kq 



q — A;p + mi 



(35) 



ME^ 



q 



2 2 , Vk^ f ^1 
m^-m + — r — 

\m^ 



(36) 



The limits of the k integral are fixed by taking cos^^ = ±1 in eq.(34) 

2 2 Vkfr 

mi — m,, — 

U — 

r^m.n.T. 



M2 



-{E^ - kr, 



m. 



2{E^ - q) 
rjk^ 



(37) 



h ■ 



M2 



(K - kr, 



2{E^ + q) 



(38) 



we solve these polynomial equations for k^ax and kmin numerically to obtain 
the kinematically allowed limits of neutrino momentum. Using these limits 
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to integrate over the neutrino momentum k we get the decay rate for pion. 
The effect of the superluminal neutrinos here is to restrict the phase space 
by restricting kmax and k^in- As a result the ratio of the pion decay width 
to the Standard Model prediction 

ro(7r ^ H = 1^ (i 

r 

is found to be — = 0.98 for = 20 GeV. 
To 

However for a 100 GeV pion the reduction can be as large as 73%. 



m 



m 



(39) 



3.2 Energy independent neutrino velocity 

We now assume the dispersion relation E"^ = (1 + 'r]')p^ in the lab frame 
and calculate the pion decay width. The amplitude squared for the process 
t^~{q) /^"(P)i^/x(fc) is 



IMP = 2GI/X 



mi-ml + v'k'('^^+2 



(40) 



The decay width is then given by 

Glf^ml f d'p d'k 



(41) 



Using the same procedure as in the last section we find the limits for the 
k integral to be 

'""^ " 2{E^ - q) ^ ^ 

" 2{E^ + q) ^ "^^ 

Solving these equations gives the following expressions for k^ax and kmm 

kmax = (E,-q+'^- A_] (44) 



kmin = V'~' (^E^ + q+'^- 



(45) 
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where 



A, 



A. 



T]'{ml - ml) 



(46) 
(47) 



Integrating over these hmits gives the decay width for pion in flight as 



ri'-' {ml - ml) (A+ - A_ - 2q) 



mf 



(48) 



In this case the reduction in pion decay width compared to the Standard 
Model prediction for an incident pion of energy 20 GeV is found to be 32% 
while for 100 GeV incident energy it can be as large as 96%. 



4 Conclusions 

We calculate the decay width forbidden process u — ?■ z/e~e^ which is allowed 
if the neutrino has a dispersion relation E"^ = m^ + + rj'p^ + {ri/M'^)p'^ 
in the context of the superluminal neutrinos observed at OPERA. We find 
that when the dispersion relation is dominated by the rj term (and the neu- 
trino velocity {vi, — 1) oc Ef,) then the mean decay length for this process 
is larger than the OPERA neutrino flight path. When the path length of 
the neutrinos is much larger than the calculated decay length then we can 
use the relation dE/dx = TE to calculate the energy loss rate. This rela- 
tion is only valid if there are multiple decays from a single neutrino over 
the path length. However for Opera neutrino energies we find that for the 
T) dominant dispersion relation the decay length is 17000 km(much larger 
than the Opera path length). In this case the energy of the neutrino beam 
< E{L) >= i?o exp(— L/r), which results in only 4 percent reduction in en- 
ergy. We also compute the pion decay width for these outgoing neutrinos 
and find that the deviation from the standard Lorentz conserving case is 2%. 

We also compute these processes assuming that the rj' term in the disper- 
sion relation dominates (and v,y — 1 is independent of neutrino energy) and 
find that the decay length for the u — )■ i>e~^e~is smaller than the OPERA 
neutrino path and this possibility can be ruled out by Cohen and Glashow 
[2]. We also find that the pion decay width is reduced by 32% and this 



9 



possibility for the dispersion relations can be ruled out as pointed out in 
[17, 22]. We have worked in framework of explicit Lorentz violation in the 
Lagrangian which gives a frame dependent dispersion relation. There exists 
other possibilities for generalising Lorentz transformation such that the mod- 
ified dispersion relations are also covariant and these theories too can evade 
the constraints of Cerenkov processes [21]. A numerical calculation of these 
processes using generalised dispersion relations has been performed in [22]. 
However in [22] the calculation was done in the centre of mass frame of the 
outgoing particles. But this calculation should be performed in the lab frame 
where eq.(3) is valid since the dispersion relation is not frame independent. 

We conclude that the OPERA neutrino measurement is only compatible 
with the energy dependent neutrino velocity, and this possibility can be tested 
at the LHC with neutrinos produced with energies above 500 GeV. 
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